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A  pseudo  two-dimensional  mathematical  model  is  developed  for  a  lithium  ion  battery,  integrating  the 
elementary  reaction  based  solid-electrolyte  interphase  (SEI)  growth  model  with  multiple  transport 
processes.  The  model  is  validated  using  the  experimental  data.  Simulation  results  indicate  that  the 
operating  temperature  has  great  effect  on  the  SEI  layer  generation  and  growth.  Under  different  charging 
-discharging  rates,  it  is  found  that  high  charging-discharging  rate  can  intensify  the  battery  capacity 
fading  process.  Different  cooling  conditions  are  then  applied  and  show  that  enhanced  surface  convective 
cooling  condition  can  effectively  slow  down  the  battery  capacity  fading.  After  that,  the  effect  of  elec¬ 
trolyte  salt  concentration  and  exchange  current  density  are  studied.  It  is  found  that  raising  the  electrolyte 
salt  concentration  can  improve  the  diffusion  property  of  lithium  ions,  and  stabilize  the  battery  perfor¬ 
mance  under  lithium  ion  consumption  induced  resistance  rising.  It  also  suggests  that  improving  ex¬ 
change  current  density  could  greatly  decrease  the  lithium  ion  battery  capacity  fading. 

Published  by  Elsevier  B.V. 


1.  Introduction 

The  rechargeable  lithium  ion  battery  is  one  of  the  most  popular 
energy  storage  devices  for  its  high  energy  density,  low  energy- 
weight  ratio,  small  self-discharging  and  environmental  suitability 
[1],  These  beneficial  properties  as  well  as  low  costs  have  established 
the  rechargeable  lithium  ion  battery  as  a  promising  candidate  for 
the  next  generation  automotive  and  aerospace  applications  [2], 
However,  during  the  charging— discharging  cycles  or  the  storage  at 
elevated  temperature  conditions,  lithium  ion  batteries  typically 
exhibit  capacity  fading  process. 

Numerous  capacity  fading  mechanisms  of  lithium  ion  battery 
have  been  proposed,  including  electrolyte  decomposition,  active 
material  degradation,  surface  film  formation,  overcharging,  self¬ 
discharging,  and  phase  change.  It  is  believed  that  the  irreversible 


*  Corresponding  author.  Tel:  +1  414  229  5639;  fax:  +1  414  229  6958. 
E-mail  address:  cyuan@uwn.edu  (C.  Yuan). 

0378-7753/$  -  see  front  matter  Published  by  Elsevier  B.V. 
http://dx.doi.Org/10.1016/j.jpowsour.2013.09.059 


capacity  loss  during  the  first  few  charge— discharge  cycles  is  pri¬ 
marily  due  to  the  formation  of  a  passivating  film  on  the  negative 
electrode  surface  [3,4],  This  so-called  solid  electrolyte  interphase 
layer  (SEI)  is  formed  at  the  first  time  when  the  electrode  contracts  the 
electrolyte  solvent  [5],  It  can  isolate  the  negative  electrode  from  the 
electrolyte  solvent  and  minimize  further  reduction  reactions  of 
electrolyte  components,  through  which  the  lithium  ions  are  trans¬ 
ported  by  vacancies  and  interstitials  to  reach  the  active  electrode. 
Typically,  the  SEI  consists  of  two  distinct  layers:  an  inorganic  layer 
and  an  organic  layer.  The  inorganic  layer  is  a  compact  layer.  It  is  close 
to  the  electrode  and  in  several  nanometers  thick.  The  organic  layer  is 
a  porous  layer  and  on  the  electrolyte  side,  which  is  usually  on  the 
order  of  hundred  nanometers  thick  [6],  The  SEI  formation  can  act  as  a 
barrier  to  prevent  the  electrolyte  solvent  cointercalation,  but  the 
excessive  SEI  formation  may  lead  to  undesirable  lithium  ions  con¬ 
sumption  and  active  electrode  area  reducing,  which  finally  causes 
irreversible  capacity  fading.  Experimental  studies  [7]  indicated  that 
the  selection  of  electrolyte  solvent  and  electrode  structure  could 


Y.  Xie  et  al.  /  Journal  of  Power  Sources  248  (2014)  172-179 


173 


have  crucial  effect  on  the  formation  and  composition  of  SEI  layer. 
Besides,  during  the  charging— discharging  cycles,  the  operating 
conditions  of  lithium  ion  battery  also  have  been  proved  to  play  a  very 
important  role  on  the  SEI  layer  growth  and  repair  [8-10],  Due  to  the 
experimental  works  cannot  fundamentally  understand  these  multi¬ 
physics  processes,  modeling  approach,  as  an  important  comple¬ 
mentary,  may  provide  significant  advantages  to  reduce  the  experi¬ 
mental  efforts  through  exploring  the  battery  performance  under 
different  operating  conditions. 

In  literature,  various  mathematical  models  have  been  proposed 
to  reveal  the  formation  and  growth  of  the  SEI  layer  [11—17],  Based  on 
the  assumption  that  the  diffusion  of  electrons  through  the  SEI  layer 
is  the  rate-determining  step,  the  parabolic  law  models  were  devel¬ 
oped  by  Peled  [11],  Broussely  et  al.  [12]  and  Ramadass  et  al.  [13], 
where  the  side  reactions  were  assumed  to  occur  on  the  interface 
between  SEI  layer  and  electrolyte.  Instead  of  assuming  electron 
diffusion  is  rate-determining,  Ploehn  et  al.  [14]  proposed  a  SEI  for¬ 
mation  model  by  considering  the  diffusion  of  solvent  through  the  SEI 
porous  layer  as  the  rate-determining  step,  of  which  the  numerical 
results  match  with  the  experimental  observations  fairly  well.  To 
study  the  detailed  surface  chemistry  of  SEI  layer  formation, 
elementary  reaction  mechanism  was  proposed  by  Aurbach  et  al. 
[15],  and  the  reaction  based  SEI  growth  models  were  developed  by 
Christensen  et  al.  [16]  and  Colclasure  et  al.  [17],  These  models  all 
focused  on  the  SEI  layer  formation  process  while  the  battery  level 
models  are  simplified.  Since  the  lithium  ion  battery  performance 
highly  depends  on  the  operating  conditions,  it  is  needed  to  couple 
the  SEI  growth  model  with  the  battery  level  model,  so  that  it  can 
provide  more  accurate  prediction  for  the  battery  cycling  behavior 
and  the  associated  processes.  However,  few  such  kind  multi-physics 
models  [18—22]  recently  have  been  presented.  Cai  et  al.  [18]  incor¬ 
porated  the  energy  conservation  into  a  multi-physics  lithium  ion 
battery  life  model  without  the  effect  of  SEI  growth.  Vazquez-Arenas 
et  al.  [19]  combined  the  SEI  induced  capacity  fading  model  with 
lithium  ion  battery  model  and  studied  thermal  effects  on  the  cycling 
process.  Awarke  et  al.  [20]  also  proposed  a  multi-physics  model  by 
coupling  the  electron-diffusion  determined  SEI  growth  model  with 
lithium  pouch  cell  model.  However,  none  of  them  has  considered  the 
complicate  reaction  mechanisms  during  SEI  layer  formation.  Since 
the  detailed  chemistry  of  SEI  formation  is  critical  to  study  lithium 


ions  intercalation  and  transport  in  the  electrolyte  solution,  the 
objective  of  this  paper  is  to  develop  a  multi-physics  lithium  ion 
battery  fading  model  by  coupling  the  elementary  reaction  kinetics  of 
SEI  layer  formation. 

In  this  study,  a  pseudo  two  dimensional  mathematical  model  of 
lithium  ion  battery  is  developed,  which  has  bridged  multi-physical 
transport  processes  with  elementary  reaction  based  SEI  growth 
model  in  a  coherent  way.  Comprehensive  simulation  studies  are 
performed  to  investigate  the  effect  of  various  operating  conditions 
on  the  battery  output  performance  and  capacity  fading  process. 

2.  Mathematical  model 

A  lithium  ion  coin  cell  has  been  utilized  as  the  physical  base  of 
the  developed  mathematical  model  in  this  paper.  Fig.  1  shows  the 
schematic  illustration  of  its  structure,  which  consists  of  a  negative 
electrode  (LixCs),  a  separator,  a  positive  electrode  (LiyM^CU)  and 
two  current  collectors  on  the  electrode  ends.  The  electrolyte  solu¬ 
tion  is  1  M  LiPF6  in  a  mixture  of  ethylene  carbonate  (EC)  and 
dimethyl  carbonate  (DMC)  with  the  volume  fraction  ratio  of  1 :2.  The 
proposed  reaction  mechanisms  of  passive  film  formation  on  the 
negative  electrode  have  indicated  that  the  inorganic  products, 
lithium  carbonate  (Li2C03),  is  one  of  the  major  passive  film 
composition  [15],  Although  the  experimental  studies  showed  that 
SEI  layer  could  also  be  formed  on  positive  electrode  [16,23],  our 
following  mathematical  model  would  only  consider  its  effect  on 
graphite  negative  electrode.  The  governing  equations  of  the  model 
include  elementary  reactions  and  multi-physics  transport  processes. 

2.1.  Elementary  reaction  kinetics 

According  to  Aurbach  [15]  and  Newman  [16],  the  elementary 
reaction  steps  of  Li2C03  formation  can  be  simplified  as  following, 

1)  The  adsorption— desorption  reactions: 


EC  +  S-EC(s) 

C2H4  +  2S-C2H4(s) 

2)  The  film  formation  reaction: 


(rl) 

(r2) 

(r3) 


CO^s)  +  2Li('s)  -  Li2C03  +  3S  (r4) 

3)  The  charge  transfer  reactions: 

Eqs)-EC(s)+e-  (r5) 

C2H4(S)  +  CO|(s)^EC(s)  +  e  +  2S  (r6) 

where,  the  subscript  (s)  stands  for  the  adsorbed  species,  S  is  the 
available  surface  site.  The  forward  and  backward  reaction  rate  co¬ 
efficients,  kf  and  iq,,  are  listed  in  Table  1.  The  net  reaction  rate  r„  can 
be  determined  as, 

r«“*?n[«f*r-*sn[‘fdr  (1) 

and  the  elementary  kinetic  description  of  charge  transfer  can  be 
formulated  by  Butler— Volmer  expression, 
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Table  1 

Elementary  reaction  rates. 

Reaction  rate  fc/(m3  mol"1  s_I)  kb  (m3  mol"1  s_1) 

rl  10s  9  x  1014 

r2  1013  10s 

r3  10s  10s 

r4  1  028  1  0-17 

r5  10”  10” 

r6  10s  10s 


Js  =  SjZf{  Ij  n  cfexp  (Z“f  Us)  -  /4  n  cj^exp  (  -  Z“f  } 

(2) 

Vs  =  <h  ~  <t>2  ~  Uf  ~  7r-Rf,  J  =  "eg,  pos  (3) 

Jneg 

where,  cfct  and  c?rd  are  the  concentration  of  reactant  and  product 
respectively;  s,-  is  the  stoichiometric  coefficient;  S,  is  the  specific 
interfacial  area;  l/Jef  is  the  reference  open  circuit  potential;  Rj is  the 
film  resistance  which  is  composed  of  initially  formed  SEI  layer 
resistance  (R/,i„i)  and  produced  film  resistance  during  charging¬ 
discharging  cycles  (R/(t)), 


Rf  =  Rf ,  ini  +  Rf(t) 

(4) 

Rf(t)  =  —■ 

Up 

(5) 

According  to  Ploehn  et  al.  [14],  the  diffusion  rate-determining 
passive  film  thickness  L(t)  can  be  analytically  given  by, 

L(t)  =  2^D0sexp(^-^jt  (6) 

here,  D(?  is  the  Arrhenius  diffusion  constant;  Ea  is  the  activation 

energy  of  the  diffusion  process;  A  is  implicit  solution  of 

,  _  ftf  m 

•Mu.co.  erfrjj  lJ 


2.2.  Transport  processes 

In  addition  to  elementary  reactions,  lithium  ion  battery  also 
includes  very  complicated  transport  processes,  such  as  mass 
diffusion  process,  charge  transport  through  both  solid  material  and 
liquid  electrolyte,  as  well  as  heat  transfer.  The  mass  and  charge 
transport  processes  are  coherently  coupled  together  at  the  active 
sites  of  electrodes.  To  elucidate  these  complicated  interactions,  the 
multi-transport  processes  are  further  modeled. 

2.2.1.  Charge  balance 

The  charge  balance  includes  the  charge  transport  processes  in 
solid  phase  and  liquid  phase,  which  can  be  formulated  using 
generic  Ohm’s  law, 

Solid  phase: 


Liquid  phase: 


a(02j-0ij)  9  /  2RT(1  - f° )  9  /  eff9(lnCj)\ 

at  oxV  jdxJ  f  a x\i  ax  J 

=  SjFJj 

0) 

where,  Cy  and  C2j  are  local  capacitances;  ipij  and  faj  are  the  po¬ 
tentials  in  solid  phase  and  liquid  phase  respectively;  t°  is  the 
transference  number  of  lithium  ion;  and  /cjff  are  the  effective 
conductivities  and  given  by  Ref.  [18], 

of  =  ffj(l  -  Sj  -  <Pj)  (10) 

K?ff  =  10-4xCj(- 10.5  +  0.668  x10-3C,-  +  0.494x10-6^ 

+0.074T  - 1 .78  x  10-5  CjT—  8.86  x  10~10 
C?T— 6.96  x  10"5T2+2.80x  10-8C/T2)2eJbruggj 

(ID 

where,  ej  is  the  porosity  and  <pj  stands  for  the  volume  fraction. 

The  local  volumetric  exchange  current  can  be  defined  as  the 
summation  of  Li/Li+  intercalation  process  (Jp  and  film  formation 
process  (Js), 

Jj=J'j+Js,neg  (12) 

Jj  =  {exp(¥^)  ■  exp(-J%) }  (13) 

Vj  =  <h-<i>2  (14) 

Where,  Uj  is  the  temperature  dependent  electrode  open  circuit 
potential,  and  can  be  approximated  by  Refs.  [18,24], 

Uj  =  Uf+(T-Tref)^i  (15) 

and  the  exchange  current  density  i°  is  given  by, 

^  =  kj(cpr-cfj^f  (16) 

In  the  modeling,  the  potential  of  solid  phase  on  the  negative 

electrode  is  set  to  be  zero,  while  on  the  surface  of  positive  electrode, 
the  current  density  equals  to  the  applied  value  /o, 

=  Io  (17) 

Since  only  electrons  are  transported  through  the  outside  load 
circle,  the  charge  flux  of  Li  ions  on  two  electrode  ends  is  zero, 


_Keff002j| 

J  9*  lx=0, 


d(h  j  -  02j) 
at 


2.2.2.  Mass  balance 

The  mass  balance  of  lithium  ions  contains  two  parts:  the  mass 
transfer  in  solid  material  and  in  the  liquid  electrolyte.  In  the  solid 
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phase,  the  mass  transfer  can  be  formulated  by  the  Fields  second  law 
based  on  a  spherical  particle  model, 


(19) 


where,  c?  is  the  lithium  ion  concentration,  D?  is  its  diffusion  coef¬ 
ficient  in  solid  materials.  Due  to  the  electrochemical  reaction  taking 
place  on  the  particle  surface,  the  consumption/generation  of  Li  ions 
can  be  formulated  by, 

ad  I 

-^er|  (20) 


But  at  the  particle  center,  the  diffusion  flux  is  zero, 

ad  I 

-D/ aF  =0  (21) 

lr=0 


Table  2 

Operating  parameters. 


Positive  electrode  thickness 
Separator  thickness 
Negative  electrode  thickness 

Faraday’s  constant 
Anode  particle  radius 
Cathode  particle  radius 
Li-difFusivity  in  anode  solid  phase 
Li-diffusivity  in  cathode  solid  phase 
Bruggeman  coefficient 
Cationic  transport  number 
Volume  fraction  of  anode  solid  phase 
Volume  fraction  of  cathode  solid  phase 
Applied  1C  charge-discharge  current 
Initial  electrolyte  salt  concentration 
Anode  solid  phase  conductivity 
Cathode  solid  phase  conductivity 
Charging/discharging  duration 
Anode  local  capacitance 
Cathode  local  capacitance 


J^moHK-1 


In  the  liquid  electrolyte,  the  material  balance  of  Li  i 
by. 


dcj  a 

ej at  =  die 


where,  D?ff  is  the  effective  diffusion  coefficient  [24], 


Assuming  the  electrode  surface  is  exposed  to  the  ambient  air, 
the  convective  heat  transfer  process  is  applied  on  the  electrode 
(22)  surface, 

-A^l  =  h(T0  -  T)  (30) 

dXlx=0,f 


£>eff  =  ig— 8.43— (54/(T— 229— 0.005c,))— 0.00022c,e( 


(23) 


3.  Solution  algorithm  and  model  validation 


The  mass  transport  for  other  species  i  is  given  by, 

‘nr  =  s(DM5!)+S>r'  l24) 

Since  the  mass  flux  is  zero  at  the  two  end  surfaces  of  battery,  the 
boundary  conditions  are  applied  as, 


-Dff9C,| 

1  9x  lx=0,l 


=  0 


(25) 


2.2.3.  Energy  balance 

The  energy  balance  in  lithium  ion  battery  can  be  formulated  by, 

pCp^r  =  ^^2  +  Qrxn  +  Qrev  +  Qohm  (26) 

where,  Cp  is  the  specific  heat  capacity;  A  is  the  heat  conductivity; 
Qrxn  is  the  total  reaction  heat  generation  rate,  Qrev  is  the  total 
reversible  heat  generation  rate,  Qohm  is  the  total  ohmic  heat  gen¬ 
eration  rate,  which  are  defined  by  Refs.  [18,24], 

Qrxn  =  SjFJj  (Vu  -  <hj  ~  Uj  -  JLRf)  (27) 

dll,- 

Qrev  =  SjFJjT (28) 


The  mathematical  model  is  solved  using  COMSOL  MULTI¬ 
PHYSICS  V4.3.  Two  submodel  geometries  are  applied:  a  one¬ 
dimensional  lithium  ion  battery  model  and  a  two-dimensional 
electrode  solid  phase  model.  As  shown  in  Fig.  1,  the  ID  battery 
model  consists  of  a  positive  electrode,  a  separator  and  a  negative 
electrode,  while  the  2D  solid  phase  model  contains  two  squares  to 
stand  for  the  solid  phase  in  two  electrodes  respectively.  The  vertical 
coordinate  of  2D  model  stands  for  the  radial  direction  of  a  spherical 
particle  while  on  the  horizontal  direction  the  lithium  ions  diffusion 
is  assumed  to  be  negligible.  The  concentration  of  lithium  ions  is 
obtained  in  the  2D  solid  phase  model,  and  coherently  coupled  with 
the  multi-physics  transport  and  reaction  processes  in  the  ID  bat¬ 
tery  model.  The  “extrusion  coupling  variables"  of  COMSOL  are 
utilized  [18], 

The  developed  model  is  first  calibrated  using  battery  surface 
temperature  curve.  The  operating  parameters  are  listed  in  Table  2. 
As  shown  in  Fig.  2,  when  the  surface  convective  heat  transfer  co¬ 
efficient  is  set  at  23.5  W  m  2  K-1,  the  numerical  results  match  with 
the  experimental  data  reasonably  well.  Then,  we  apply 
23.5  W  m-2  K-1  on  the  battery  surface  and  simulate  the  battery 
discharging  process  on  the  first  cycle,  the  model  predictions  also 
match  with  the  experimental  data  fairly  well.  The  developed 
mathematical  model  thus  can  be  used  for  further  numerical 
studies. 

4.  Results  and  discussion 

4.1.  SEI  layer  growth 


n  rr  f9' M54-,  (d^>2\2  2k^KT7  Q \  9(lnc) 6</>2 

Qohm  =  (HarJ  +Keff^J  +—^{1-t+)-aiTW 

(29) 


The  SEI  layer  generation  process  is  greatly  affected  by  the 
applied  electrode  material  and  electrolyte  solvent,  as  well  as  the 
variation  of  battery  operating  conditions  [16].  Since  the  compo¬ 
sition  of  electrodes  and  electrolyte  solvent  are  definitely  known, 
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Fig.  2.  Comparison  between  numerical  results  and  experimental  data. 


the  effect  of  different  operating  conditions  become  more  impor¬ 
tant,  one  of  which  is  the  operating  temperature.  Fig.  3  has  simu¬ 
lated  its  effect  on  the  SEI  layer  growth.  Three  initial  temperature 
conditions  of  15  °C,  25  °C  and  60  °C  are  applied  and  the  battery 
discharge/charge  rate  is  set  at  1C.  Numerical  results  show  that  the 
thickness  of  SEI  film  keeps  increasing  when  battery  is  operating, 
and  a  higher  initial  operating  temperature  clearly  leads  to  the 
generation  of  a  thicker  SEI  film.  Compared  with  the  electro¬ 
chemical  analysis  of  Christensen  et  al.  [16],  the  predicted  SEI  film 
growth  is  a  little  slower.  This  may  be  because  our  model  has 
considered  the  effect  of  multi  transport  processes  in  lithium  ion 
battery  as  well  as  the  electrode  microstructures.  Since  SEI  layer 
growth  can  cause  undesirable  lithium  ions  consumption  and 
resistance  rising,  in  the  following  sections  we  will  study  this  effect 
on  the  global  cell  performance  under  different  operating 
conditions. 


4.2.  Effect  of  charging-discharging  rate 

The  charging-discharging  rate  has  significant  effects  on  the 
multi-physics  transport  processes  of  lithium  ion  battery.  Three 
charging— discharging  rates  of  1C,  0.75C  and  0.5C  are  studied  in 
Fig.  4.  The  initial  operating  temperature  is  set  at  25  °C.  As  shown  in 
Fig.  4(a),  the  cell  capacity  of  0.5C  case  is  much  higher  than  the  other 
two  cases  at  the  same  cell  voltage  on  the  first  cycle.  It  can  hold 
around  14  Ah  m-2  at  3.7  V  while  that  of  1C  case  can  only  achieve 
7  Ah  m-2.  This  is  due  to  high  charging— discharging  current  density 
can  intensify  the  electrochemical  reactions  as  well  as  the  SEI  layer 
generation  reaction,  and  lead  to  the  rising  of  cell  polarization 
resistance.  Fig.  4(b)  illustrates  their  effects  on  the  battery  cycling 
performance,  where  both  the  charging  and  discharging  period  are 
set  to  be  3000  s.  Upon  three  current  conditions,  the  voltage  drops 
during  discharging  process  are  all  kept  increasing.  Moreover,  it 
shows  that  the  fading  trend  of  1C  case  is  much  worse  than  the  other 
two.  This  observation  indicates  that  high  charging— discharging 
rate  can  greatly  intensify  the  SEI  growth  induced  capacity  fading 
process  during  battery  cycle  operating. 


4.3.  Effect  of  cooling  condition  and  initial  operating  temperature 

The  effect  of  different  cooling  conditions  on  the  cell  temperature 
variation  is  shown  in  Fig.  5.  The  discharge  current  rate  is  1C  and  the 
initial  operating  temperature  is  set  at  25  °C.  Obviously,  the  cell 
surface  temperature  keeps  rising  during  the  discharging  process. 
When  the  surface  convective  heat  transfer  coefficient  is  increased 
from  zero  (insulation  boundary)  to  50  W  m  2  K,  the  cell  tempera¬ 
ture  ascending  trend  is  clearly  restrained.  It  is  decreased  from  as 
high  as  75  °C  to  around  35  °C  at  10  Ah  m-2. 

Since  cooling  condition  induced  temperature  variation  is  sig¬ 
nificant,  Fig.  6  has  illustrated  its  effect  on  the  battery  capacity 
fading,  where  the  discharging/charging  rate  is  set  at  1C  and  the 
discharging  and  charging  period  are  both  3000  s.  Two  convective 
heat  transfer  boundary  conditions  of  h  =  0  W  m-2  K  (insulation) 
and  h  =  23.5  W  m-2  K  are  compared  in  Fig.  6(a).  As  shown,  the 
battery  capacities  of  both  two  cases  keep  reducing  from  the  1st 
cycle  to  9th  cycle,  but  it  is  obvious  that  the  capacity  fading  of 
insulation  boundary  case  is  much  faster  than  the  better  cooling  one. 
It  means  that  well  convective  cooling  condition  can  effectively  slow 
down  the  lithium  ion  battery  capacity  fading  process.  Similarly,  the 
effect  of  initial  operating  temperature  on  the  capacity  fading  at  9th 
cycle  is  studied  in  Fig.  6(b).  It  shows  the  higher  operating  temper¬ 
ature  is  applied,  the  faster  fading  process  will  be  generated. 


4.4.  Effect  of  electrolyte  salt  concentration 

The  electrolyte  salt  concentration  is  another  important  factor 
determining  lithium  ion  battery  performance.  In  Fig.  7,  the  con¬ 
centration  profiles  of  the  binary  electrolyte  salt  in  discharging 
process  at  different  cycles  are  simulated.  The  initial  operating 
temperature  is  set  at  25  °C  and  the  discharge  current  rate  is  1C.  It 
shows  the  electrolyte  concentration  profile  of  6th  cycle  is  much 
flatter  than  that  of  other  cycles,  which  suggests  that  the  diffusion 
property  of  electrolyte  solvent  on  6th  cycle  is  better.  This  is  because 
more  heat  will  be  generated  as  the  charge-discharge  cycle  goes  on, 
which  can  benefit  the  diffusion  process.  Besides,  it  can  be  noticed 
that  the  concentration  profile  patterns  of  6th  cycle  and  4th  cycle  are 
very  close.  Because  the  SEI  layer  thickness  keeps  increasing  during 
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the  cell  operating,  this  observation  suggests  its  passive  effect  on  the 
electrolyte  solvent  diffusion/reaction  is  also  increasing. 

According  to  the  elementary  reaction  mechanism,  the  step  re¬ 
action  rates  are  determined  not  only  by  the  applied  battery  oper¬ 
ating  conditions  but  also  by  the  concentration  of  each  reactant. 
Fig.  8  illustrates  the  effect  of  different  initial  electrolyte  salt  con¬ 
centration  on  the  capacity  fading  of  lithium  ion  battery,  where 
three  concentrations  of  0.001  mol  cm-3,  0.002  mol  cm-3  and 
0.004  mol  cm-3  are  employed  respectively.  The  initial  temperature 
is  25  °C  and  other  operating  conditions  are  the  same  as  above  cases. 
The  battery  discharging  performances  are  compared  at  the  4th 
cycle.  It  can  be  seen  that  the  battery  with  0.004  mol  cm-3  elec¬ 
trolyte  salt  has  the  slowest  fading  process.  This  can  be  explained  by 
that  adding  more  lithium  ions  to  the  battery  electrolyte  solvent  can 
improve  the  diffusion  property,  and  stabilize  the  cell  performance 
under  lithium  ion  consumption  induced  resistance  rising.  Besides, 
different  cooling  methods  on  three  electrolyte  cases  are  also 
compared  in  Fig.  8,  which  shows  that  enhanced  cooling  condition 
can  obviously  slow  down  their  capacity  fading  process. 

4.5.  Effect  of  exchange  current  density 

Theoretically,  the  exchange  current  is  the  ongoing  current  of  a 
redox  reaction  in  both  cathodic  and  anodic  directions,  which 


Fig.  6.  Effect  of  a)  different  cooling  conditions  and  b)  initial  operating  temperature 
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Fig.  7.  The  concentration  profile  of  the  binary  electrolyte  salt  in  different  cycles  during 
discharging  process. 


reflects  the  intrinsic  rates  of  electrochemical  processes.  It  depends 
critically  on  the  nature  of  the  electrodes  and  electrolyte  solvent 
including  the  physical  structure,  chemical  composition  and  the 
electrochemical  properties.  Since  it  has  great  influence  on  the  cell 
performance,  Fig.  9  has  illustrated  different  exchange  current 
conditions  during  the  battery  cycle  operating.  To  make  a  compar¬ 
ison,  the  exchange  current  densities  are  chosen  to  be  50%  and  5 
times  of  its  base  value  respectively.  The  cycle  process  is  simulated  at 
25  °C  and  1C  charge/discharge  rate.  It  is  clear  that  the  lowest  ex¬ 
change  current  density  case  (0.5  Io)  has  the  fastest  fading  process. 
Since  the  exchange  current  density  reflects  the  electrochemical 
reaction  rate,  it  can  be  concluded  that  the  method  like  increasing 
the  active  surface  area,  raising  Li  ion  concentration  that  can 


process. 


intensify  the  electrochemical  reactions  could  have  positive  effect  to 
slow  down  the  SEI  growth  induced  battery  fading  process. 

5.  Conclusion 

A  pseudo  two-dimensional  mathematical  model  of  lithium  ion 
battery  is  developed  by  coherently  coupling  the  multi-physics 
transport  processes  with  elementary  reaction  based  SEI  growth 
model.  The  model  is  validated  with  experimental  data.  Simulation 
results  indicate  that  the  operating  condition  of  lithium  ion  battery 
has  significant  influence  on  the  SEI  layer  generation.  Under 
different  charging-discharging  rate,  the  comparison  study  shows 
that  high  charging-discharging  rate  can  greatly  intensify  the  bat¬ 
tery  capacity  fading  process.  With  changing  surface  cooling  con¬ 
dition,  the  numerical  results  also  show  that  well  convective  cooling 
condition  can  effectively  slow  down  the  battery  capacity  fading. 
The  battery  performances  of  three  different  electrolyte  salt  con¬ 
centrations  are  then  simulated,  which  suggests  raising  the  elec¬ 
trolyte  salt  concentration  can  improve  the  diffusion  property  of 
lithium  ions,  and  stabilize  the  cell  performance  under  lithium  ion 
consumption  induced  resistance  rising.  The  effect  of  exchange 
current  is  finally  simulated.  It  is  found  that  the  method  of  raising 
exchange  current  like  increasing  the  active  surface  area  or  raising  Li 
ion  concentration  could  slow  down  the  battery  fading  process 
through  intensifying  the  electrochemical  reactions. 
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List  of  symbols 


cfrt,  cfrd  the  concentration  of  reactant  and  product 

Cp  the  specific  heat  capacity 

D|ff  the  effective  diffusion  coefficient 

0(3  Arrhenius  diffusion  constant 

£fl  the  activation  energy  of  the  diffusion  process 

F  the  Faraday’s  constant 
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i?  the  exchange  current  density 

/o  the  applied  charging— discharging  current  density 

Jj  the  local  volumetric  exchange  current 

kf,  kb  forward  and  backward  reaction  rate  coefficients 

Rf  the  film  resistance 

Si  stoichiometric  coefficient 

Sj  the  specific  interfacial  area 

tfj.  the  transference  number  of  lithium  ion 

U]  the  electrode-open  circuit  potential 

l/Jef  the  reference  open  circuit  potential 

0ij,  <hj  the  potentials  in  solid  phase  and  liquid  phase 

of,  Kjff  the  effective  conductivities 

A  the  heat  conductivity 
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